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This presentation will describe the effort at Sandia National 
Laboratories Albuquerque with emphasis on the methodologies and 
techniques being used to generate strict hexahedral finite element 
meshes from a solid model. We utilize the functionality of the modeler 
to decompose the solid into a set of non~ intersecting meshable finite 
element primitives. The description of the decomposition is exported, 
via a Boundary Representation format, to the meshing program which 
uses the information for complete finite element model specification. 
Particular features of the program will be discussed in some detail 
along with future plans for development which includes automation of 
the decomposition using artificial intelligence techniques. 


★This work performed at Sandia National Laboratories was supported by 
the U.S. Department of Energy under contract DE-AC04-76DP00789. 
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A proposed system to generate finite element models directly from 
the solids model database Is presented. This system includes *uto- 
matic error analysis with adaptive gridding for equilibration of the 
error estimator in use. The complete specification the finite 
element model including boundary conditions and material identifiers 
is produced to a neutral output file. An illustrative example depict- 
ing*^^ the state of implementation of the proposed system is contained 
within. Current research is also briefly described. 

Introduction 

The advancing technology of computing hardware and software is well 
represented by the current Computer Aided Design (CAD) systems employ- 
ing solids modeling. These solids modeling systems, u «der development 
by both universities and industry, have the obvious benefit for the 
realistic visualization of three-dimensional (3-D) objects. The mos 
important benefits of solids modeling, however do not lie in the 
solid model itself, but in the subsequent applications which utilize 
the valid and unambiguous geometric information available. In other 
words, the advantage of solids modeling is not as a stand alone appli- 
cation but as a means of creating a geometrical database to unify a 

number of applications. Indeed, users and vendors currently seem to 
be concentrating their efforts at integrating the solids model data- 
base in the areas of Finite Element Modeling (FEM) and Numerical 
Control (NC) Programming. Solids modeling does appear to have the 
potential for unifying the design, engineering, and manufacturing 
areas of industry. 

At Sandia National Laboratories a unified geometric database is 
expected to reduce design time and yield added reliability and 
optimization of the designed systems. A Joint effort between the 
Engineering Sciences and the Computer Aided Design Departments has 
been defined and is being pursued to integrate the Computer Aided 
Engineering (CAE) activities of the Engineering Sciences Department 
into the automated design and manufacturing process. TCie primary 
vehicle for this effort is the utilization of improved model genera- 
tion capabilities with emphasis on advanced geometric definition and 
automatic mesh generation for FEM. In particular, the utilization of 
the CAD geometrical data and hence the elimination of the error prone 
reentry of such data is considered essential. 


♦Member of Technical Staff. CAD Technology Division 2814, 
Sandia National Laboratories, Albuquerque, NM, 87185. 
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Thl« nancr describes the effort underway et Sandia for 
integration ? of FEM Mesh generation utilizing PADL-2 . the 

Constructive Solid Geometry (CSC) system produced at the University of 
Rochester. In general, because of the commercially •v-il^le end 
locally developed finite element analysis codes in use D *JJ;*«iJt - 
requirement for the use of hexahedral elements in 3-D FEM exists. 
This coupled with the large number of nonlinear finite *l®“«nt 

analyses performed, prohibits us from considering the automatic 
tetrahedralization work [CAVE85] developed at General Motors Research 
Labs or the modified-octree work [YERR84] performed at R ?"***l a f* 
Polytechnic Institute. The finite element mesh generation philosophy 
we are pursuing is divided into two primary phases, 1) initial mesh 
installation utilizing the available CAD geometric data base *nd 2) 
mesh optimization including mesh improvements based on R eo “ e * r *J a ?; 
aspects" of the initial mesh and automatic error analysis coupled with 
Sode griding techniques to obtain uniformly reliable answers 

throughout the domain of analysis. 

The following sections of this paper describe in some detail the 
relevant topics including 1) solids modeling, 2) application inter - 
face 3) initial mesh generation, 4) mesh improvements, 5) and error 
analysis and adaptive grldding. An illustrative probl.a depleting the 
state of implementation of these topics is included. 

Solids Modeling - A Geometrical Basis for Applications 

The classical geometrical CAD database is the so-called 
•wireframe- format. To define wireframe, we introduce the notion of 
an ed-e. For us, how an edge is actually represented within soue 
iom P !t« database is unimportant. Only the idea that an edge results 
from the intersection of two distinct surfaces matters. An edge is 

one -dimensional in a parametric sense. That is*. .J^oieh 

(x v z) on an edge is in Euclidean 3-space it can be derived through 
a* system of equations depending on only one independent parameter of 

the form 

x - X(s) y - Y(s) z - Z(s) . 

He-e y Y and Z are functions of the independent parameter s which is 
bounded in* the closed interval A wireframe representation 

then models a solid by simplistically specifying certain edges of the 
solid. Typically, those edges defined for a given solid correspond to 
the bounding edges of the domain being considered. "jj'j 

specification, referred to as instancing, is accomplished through a 
choice of a particular type of edge (say a line or circular *eg»ent) 
with a rigid motion and any other necessary parameters (say curvature) 
to complete the definition. 

New geometrical modeling technologies are becoming popular. The 
two most popular technologies are CSC and Boundary Representation (B- 
Ren) CSG systems define solids as Boolean operations (union, dif- 
ference .intersection) of simpler primitive solids (blocks, spheres 
wedges, cones) instanced by size and location. The B * R *P’ . 

other hand, is a heirarchical extension of the wireframe form st. In 
the B-Rep, solids are described as a collection of instanced 'ky t^P , 
size, and location) faces, each of which in turn are composed by a 
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number of «dges. Explicit mathematical descriptions °* b ° tb ****! 

•~5 .d«. Ire usually available. The user interface for the CSC and 
a Ren Sterns Ippear to be unifying with each other borrowing from the 

rcr «od*ler Is found in several B-Rep modelers. Similarly, a sweep 
HI SSiS tioi of edges to create faces and the sweeping of faces to 
form volumes have begun to appear in some CSG formulations. 

We are of course, interested in utilizing an unambiguous and 

r is? sr,r. zz 

end only one solid. We do not have strict one-to-oneness since there 
is no unique representation for a given solid, but only 

for e^ry representation. Indeed, in any of the currently available 
geometrical modeling systems, there is no unique representation for a 
■olid There are as many definitions of a solid as there are users. 

-* ... » #rm "valid" we imply that for any representation we 

Strict rliJld -Uhougt It need not b. r..Ut.bl. fro. . 

manufacturing point of view. 

It is easy to imagine that wireframe representations are neither 
unM ,blguo« *nor velld. Indeed , there ere . .yrl.d of coun or- 

4va _ nl r c testifvinE to this. On the other hand, both CSG and B Rep 
systems have the ability to produce unambiguous and valid descrip- 
tions. Our work in automatic finite element analysis h ** b ** n 
on Jhe unambiguous and valid geometric description available within 
the CSG modeler PADL-2. The choice of PADL-2 has been more a ®?tter 
of convenience , .Inc. the .cure. cod. end exp.rtl.e ere .v.llebl. « 
c.ndla than a matter of preference of CSG over B-Rep. In fact, it 
be arsued that the B-Rep facilitates certain applications, for 
example FEM and NC programming, that are primarily surface oriented. 

For our implementation of the mesh generation we use the B-Rep, 
as supplied through a conversion routine available in PADL- . 
conversion routines are generally well understood and details of the 
pISl- 2 implementation can be found in [HART81]. Our development thus 
is considered generic in the sense that any solid modeler capable of 
, . at-oiv daliverine a B-Rep independent of its own internal repre- 
:USS£ “.m! !?«»!*! the re develop- 

ing. 

We have realized the benefits of using a valid and unambiguous 
solid model as neither the geometry nor the topology h * s 5 ° b * *“ p * 
elemented For wireframe applications it is quite typical that either 
additional topology or geometry has to be supplied before applications 

are undertaken. 

Application Interface - The Link Between Geometry and Applications 

The idea for using arbitrary solid modelers in conjunction with 
various applications is known under the broader category of 
-application interface". An application interface has been likened to 
a "software bus" enabling applications to communicate directly to 
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•olid* modelers for the purpose of Interrogeting or modifying the 
solid model. To dote no etenderd epplicetion interfece exists for the 
available solids modelers although efforts [CAMI86] to this end have 
been underway for some time. Still, some solids modelers make 
available to some extent the modeling operations required by locally 
de/e loped applications. Application interfaces can be thought of as 
part of the solids modeler which make the internals of the modeler 
transparent to the application. 

Ve have been able to use a number of the available routines 
within PADL-2 to facilitate the interface to the finite element mesh 
generator. These include routines for identifying and utilizing the 
geometric entitles within the representation. For example, routines 
pertaining to the storage management structure, the rigid motion 
facility, and the computational geometry package have been used to 
discretize the body for mesh generation. Other utilities necessary 
for linking our application to PADL-2 have had to be defined and 
developed. These include routines that format the B-Rep available in 
PADL-2 for export to applications and the corresponding routines to 
read the representation into the mesh generator. The mesh generator 
also requires contiguous lists of edges and faces, called loops, which 
PADL-2 does not require. These have been developed. Redundant edges 
and faces are either necessary or add to the robustness of a aolids 
modeler, but are detrimental to mesh generation. Algorithms to 
identify and eliminate redundant faces and edges have been 
implemented. Finally, although PADL-2 contains routines usable to 
discretize edges, none existed to discretize surfaces. 

The above development has allowed the mesh generator to directly 
create finite clement models from a solids description while guaran- 
teeing that all nodes defined for the mesh either lie in the body or 
are exactly on the surface. This group of routines are necessary 
-within another solids modeler for our implementation of the finite 
element mesh generator. 

Mesh Generation - An Application for Solids Modeling 

In this section the philosophy for generating hexahedral finite 
element meshes from a solid s model database is presented. We proceed 
by briefly describing one technique for generating hexahedral meshes 
that is representative of the classical methods used. The method we 
are pursuing for mesh generation is an extension of these ideas im- 
bedded in new technologies, namely solids modeling and feature recog- 
nition. 

A hexahedral mesh can be constructed through a coordinate trans- 
formation in conjunction with higher order approximating functions. 
More specifically, the geometry of the body is constructed using 
hexahedral subregions each having six well defined faces and twelve 
edges. The description of each hexahedron requires the coordinates of 
eight corner points and one interior point of each edge for a total of 
twenty points. During the construction of a particular hexahedron, 
faces which are coincident with previously defined faces are identi- 
fied. Thus, coincident nodes on coincident faces are assigned the 
same node number. Finally, a consistent number of divisions along 


orthogonal* directions for oach haxahsdron Is 


thrsa "mutually 
specified. 

. . Ksvahedrel elements can then be Installed In each 

** , _ | / f \ ,A^ r • t<11 via naps aiven by x m A(r,i,t) , y 

S* where S “ [C * » * — • * jj ^ the usually polynomial functions X, 

obtained In that portion of the subregion. 

r S? srS-irrilnH 

^ clarification. . cap. 

^ilitiy of the solids modeler PADL-2. 

In the discussion of a classical hexahedron mesh §« ner J^‘ 
da.crlSad th. s.oj.trlc d.fl ? Jtlon of Cf th^body “ 

S' construction la* 

^U."°^: r '«on.^tr can b. overly tin. conauning for .11 but 
the simplest cases. 

*7 °xh. *& “:r\rn^ c 'i:ur.:d:irr°r. t ^iu::d lon oir 

p““« 0 Phy for cubr.glon ‘^.‘body*!™'^!" of 

?“•? ?*£ °^i ch ~-' ‘i«sir.;sr “iStS?.^ 

installed. ... del * s decomposed using the primitives and 

Boolean* * operations of ®Ua aollda ncdeler lnto^ att of f ini te^lanant 

finite elemen p toDolologically equivalent entities. For 

exiSpl^ eS any Volume defined by one surface, topologically equivalent 
to a sphere, will be able to be meshed. 
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Allowing wore general finite element primitive* either 
necessitates the definition of new mapping techniques or a decomposi- 
tion of each of the finite element primitives into a collection of 
hexahedra. This last alternative is easily accomplished. Figure 1 
shows the decomposition of the standard geometric primitives. This 
decomposition of the finite element primitives will be automatic in 
the solid modeler and transparent to the user. 

The mesh generation is only automatic in each finite element 
primitive. Presently human interaction is required for the primitive 
decomposition. Work is beginning in the area of feature recognition, 
as applied to recognizing the finite element primitives, to automate 

this process. 






Figure 1. Decomposition of the Standard Geometric 
Primitives into Hexahedron. 


Mesh Improvements - Assuring Geometrically Cood Meshes 

In the previous discussion of the mesh generator we did not 
enumerate the characteristics of a "good" mesh. We do so now. Some 
characteristics of a good mesh are 1) gradually changing element 
sizes, 2) gradually changing element shapes, and 3) as nearly rectan- 
gular (even cubical) elements as possible. These characteristics have 
Important numerical consequences. For example, the third condition 
assures us in practice of a well defined (one-to-one and onto) coor- 
dinate transformation during the stiffness matrix formulation. In 
addition, all the above characteristics attempt to maintain the condi- 
tion numbers of the stiffness matrices generated for two nearby ale- 
merits to be similar. 
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Our approach to generating good meshes i» an extension the 
ideas lncoporated in the two-dimensional (2-D) mesh generator QMESH 
[JONE74] . Only the necessary details o f these developments X^im* 
-iven in thi* section. A more complete description of the 2 
pleaentation can be found in the QMESH documentation. 

In our aesh generator, like QMESH. the initial aesh i««valuated 
and iaproved through a series of processors working in tandea. The 
processors have the capabilities to autoaatically reposition nodes, 
delete eleaents, and rearrange the topology in an attempt to improve 
the element geometry. The aesh improvements, as we now discuss, are 
only concerned with the geometrical aspects of the aesh. 
ciency of the aesh with regard to accuracy is discussed in the next 
section of this paper. The general concepts of the algorithms for 
node smoothing, topology restructuring, and element deletion are now 

described. 

The node repositioner, or smoother, consists of attempting to 
have the nodes equidistant and the elements having equal volumes. 
Requiring the nodes to be equidistant is tantamount to requiring that 
every node is at the average location of all its neighbors. Symboli- 
cally we have 


- £<IvlJVl*i> 


where n is the number of neighboring points each with coordinate (Xj, 
y z ). This formula is the one applied in the smoothing code but 
in* a slightly altered fashion. The expression is rewritten as 

(x,y,z) - + n^ x i‘ x 0 ^ *^ y i' y 0 ^ '^ Z i’ Z 0^ 

or more succinctly 

(x t y,z) " (xQ.yo^o 5 + 

Here lx v z ) snd (x,y,z) are respectively the old and updated 
positions’ of °the node being moved and V L - is the "Laplacian" 

movement vector. 

Only a related form of the volume equilibration has been 
considered to date. Instead of requiring the volumes for all elements 
to be the same, we impose that each of the areas of the faces onfall 
the elements are equal. This has been proposed to more fully utilize 
the capabilities already developed in QMESH. This requirement is 
represented in the formula 

v. -y A ir^ife v 

A if +A ib 

where V is the "Area-Pull" movement vector (corresponding to in 
the node A equi distribution) applied to the node in question. The A f 
and A v refer to the areas of the face "in front" and "in back of 
the nlde. Again, more complete description of the formulation in the 
two-dimensional setting can be found in (J0NE74] . The Laplacian and 
Area -Pull moment vector for a node are incorporated through a convex 
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combination of the two. That la the moment vector for a given node is 

taken as 


V - oV A ♦ <l*a)V L 

with ae (0,1] e user selectable parameter. 

The next capability is the restructuring of the element topology, 
i.e., the element connectivity. By this we mean the process of eras- 
ing an interface plane and drawing it differently to improve the 
geometrical shape of the neighboring elements to the plane. To assess 
the element shapes, three element evaluator functions referring to the 
angle condition, the aspect ratio, and the product of these two have 
been defined. These definitions are extensions to those developed by 
Jones for QMESH. 

The operation of the restructuring process is then the following: 
the condition numbers for all the elements are evaluated and a list of 
the twenty- five worst (largest) is saved. The processor attempts to 
Improve the worst element in the mesh. If no improvement is made in 
any of the first ten worst elements, the processor quits. If a re- 
structure is accomplished, the list of worst elements is updated and 
the process is continued so long as a restructure is performed among 
the ten worst elements of the mesh. 

The final processor contained is the element deleter. This 
processor attempts to improve the mesh by deleting elements. Element 
deletion is similar in nature to the restructuring processor. This 
processor sweeps through the mesh to make a list of the five worst 
"rhombic" elements. The measure of how rhombic an element, termed the 
R-number, is defined as the ratio of the length of the shorter diago- 
nal to the length of the longer diagonal. If the R-number of an 
element is less than tan(V/2), where V is normally forty-five degrees, 
the element is placed in the candidate list for deletion in ascending 
order. The more rhombic an element is the smaller the R-number. The 
tolerance parameter tan(V/2) is the R-number of a parallelogram with 
opposite angles of V and is not simply a measure of how sharp an 
element is. The program then starting with the worst (smallest) R- 
number, attempts to eliminate the element. As soon as an element is 
deleted, control returns to the calling program. 

The sequence that the processors operate on can effect the out- 
come of the mesh. A general method for specifying the sequence of the 
processors has been implemented. The entire sequence is iteratively 
performed until convergence (no more node smoothing, element restruc- 
turing, or element deletion) is attained. 

The full capability of the mesh improvement has not been imple- 
mented to date. Only those capabilities corresponding to each of the 
lofting planes generated in the primitives are acted on. We can show 
that for certain limited cases this is a partial implementation of the 
entire algorithm for a true 3-D mesh improvement. We have observed 
that the primary processor functioning is the smoother, attributing to 
the Initial quality of the meshes generated. 
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Error ta.ly.l. .nd Bd.pt lv. Crlddtng - Cen.r.tlng CompototlonoUj 
Opciama Methtt 

Th* tooics of error analysis end adaptive gridding have attracted 

sr. ia S 5 Si 4 "sSS; 

S* “r“ r «n *TSdSid I 5 l£*Si , I-Kr “ dU^'/ef 'freedom 

«• gjssr-- y^ss^iTSTfat?:^ 

b7 lnor.«lng th. "“»« referrtd to t h. -p-method", 

space. decree of the polynomial used over each 

Involve. ,' h * th f number of^l.ment. In th. do..in 

fixed" Thl. refinement, .. U.U.II7 Implemented, h*. no «£<••« °? **“ 
UprUl.IUU of th. geometry .nd hence th. nltl.l .e.h .hould be 
developed to Include .11 th. Impotent geometrlc.l tcpecte. 

The ..cond method of lncr...lng the .ccur.cy of th. finite ele- 
ment method is termed the -h-method". In the h-method, the degree ®‘ 
the approximating functions over each element is ^inteined while th 
Bortion of the domain spanned by each element is decreased. I 

Sords, additional elements (and thus nodal P oi " c *> *5* no im- 

. ’ A-«<n in oractical applications of the h*method, no im 

provement of the geometrical model is attained. The ^f^etica 
aspects of the h-method are probably more developed than thos 

pertaining to the p-method. 

The difficulty with both the above approaches is that extensive 
, 4 P Vl m r<m .r. reouired for the finite element analysis codes in 
SJ UdIJ™o take advantage of these developments. A * h J rd 
tive for adaptive gridding is possible. Errors in numerical methods 
Ur polntvlse^ dependent . Th.t 1., error. In .n .n.ly.1. u.u.lly v.ry 

from P one element to .nother. One r...on.bl. go.l t0 *" £v * f? r j!rv 
uniformly reliable answer with the same error associated with every 
element ^ Thus if elements are concentrated in the are. where errors 
Ire Urge while decreasing the number of elements where «rors are 
small we can hope to produce such a result. In effect, we are tryi g 
to automate the technical expertise applied by computational stress 
analysts in achieving reasonable results. 

The methodology currently being pursued is the 
orimarily because of the large investment in conventional FEM tech 
£ 5 ;® not involving the h- or p-methods of refinement. It may be 
Ur. .0 ll UmU Ufor. .oftv.r. P U re.dlly .viable for .pplylng 
eeneral h- or p-methods especially for the non-linear problems of 
interest at Sandia. The question now arises: HOW do we introduce the 
adaptive grading techniques in our calculations? To answer this we 
lolk al thf composite partis of the problem. They are error evaluation 
and node distribution. 

The problem of finite element error evaluation has been «tudied 
extensively. Sophisticated theoretical work has been done by Babuska 
.„d Rhelnboldt ( BASU78. , BABU78b , gABUSO) ; K. Mlll.r .nd «. Mlll« 
[MILLS la] ; K. Miller (MILLSlb); and Babuska and A. Miller on 

error evaluation. To paraphrase their work without extensive techni- 
cal details, an estimate of the local error in an energy norm at a 
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given node is derived by considering the surrounding elements to thet 
node. Here the energy norm for e function f Is defined es 


imi E(D , - J D i vf i *-• 

The Indicated Integration Is carried out over the domain covered by 
the neighboring elements of the node of interest and is notated >y D. 
It is assumed that the finite element solution, denoted by u^, to the 
desired partial differential equation (PDE) is available. A model 
problem, perhaps mimicking the actual PDE being solved, with boundary 
data corresponding to u. is then solved on D and denoted by w. It is 
then reasonable to assume that the quantity ll u h* w * *E(D1 V 

approximates the local energy error 1 1^ • ullp/m ^e*?e u li the 
exact (and usually unavailable) solution. Varliui refinements and 
extensions to this idea are the topics of the references cited. 


The error indicator we are currently experimenting with is 
different from the one presented above. Our error indicator is 
simpler to evaluate and attempts to estimate the maximum pointwise 
error in an element rather than a local energy error. Briefly, under 
conditions which are usually satisfied, it can be shown that in the 
maximum norm the finite element solution is the optimal solution 
available from the span of the basis functions. Thus 


||u-Uhll L . < C min ||u-X|| Li# 
XeS h 


where II * I |. denotes the maximum norm and S. denotes the finite 
•lenient eubspace depending on the choice of the discretization h and 
the approximating polynomials S used. Then, on an element k approxi- 
mation theory yields 


I Il» (K) < 


CVI 7 U|I L- (K) 

:«4i (k> . 


2 

Here D u denotes 
diameter of the 
first inequality 
second inequality 
norm is currently 


the generic second derivative while h^ is the 
element. In the case of linear elements we use the 
replacing V u by ?u. . Eor quadratic elements the 

can be used replacing D\i with D*u. . The maximum 
estimated at the quadrature points or the elements. 


This error indicator is conservative for problems with smooth 
solutions. It does not take into account the full order of the 
polynomial approximation for such problems. It is known that for 
problems with 0 smooth solutions, the error using linear elements would 
be of order h Z and not.h. Similarly for quadratic elements the error 
would be of order h J and not h . This error indicator may be more 
suitable for non* smooth problems such as in shock calculations. 
Again, the primary advantage of the proposed error indicator is its 
computability and it is a pointwise estimate. 


Ve now consider the problem of how to distribute the nodes to 
equilibrate the error indicator and obtain a uniformly reliable 
solution throughout the domain of anyalysis. This problem is easily 
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addressed with the development of the 8«o“*trlcel 
cussed earlier. As noted, the moveaent vector as 
Area- Pull and Laplaclan smoother Is given as 


smoothing as dis- 
considered In the 


V - «V A *(l-o)V L 

It is reasonable to expect that the new node movement vector defined 

as 


V - aV A e A ♦ (l-a)\« L 

. j _ rcnrtient the errors associated with the Area-Pull 

and f *the A Lap lac tan movement vector respectively, would cluster the 
5LJ5 vhera.rrors .r. Ur,.. For.xa.pl., ® A could b. d.flH.d .. 

*L COU .. i Overall the node distribution for adaptive grading 

*u r »un « of th. g.o«trlcal nod. nothin,. 

An Example - Current Capabilities 

As an illustration of the procedure for creating a finite element 

„„ h .„d S2ST *%£ 

nrooerlv instanced blocks and the difference of a properly Instanced 
cylinder ^e solid is decomposed into six finite element primitives. 
Xl Solids model of the decomposed L-br.cket is shown in Figure 2 
Si f our S p r im i t i ve s surrounding the hole were created by first 
defining .Coordinate system with respect to the axis of the hole A 
properly instanced wedge with it. apex P^all* 1 to the holes J* 1 * 

kr.pp»~«-ni^ 

instanced blocks with the L-bracket. 

The B-Rep of the decomposed bracket was transferred to the mesh 
> «»veral of the routines imbedded in the application 

Discretization data supplied via the user interface in the mesh gene 
»Iul“d in rh. ..sh Indicated In Flgur. 3. *PP<-opr .t. 
rl.l Indicators .nd boundary condition flags »«. 
within th. FEM 

sss i.d £o nnn. fo :i.::« y ‘i‘oi.i ««. s *f » w 

linear elastic static analysis was accomplished. Typical results 
shoving the minimum principal stress contours is shown in Figure 3. 


Fixure 4 shows the result, of the automatic error analysis. 
Areas of large error are indicated as a high density of error contours 
and* Corresponds to locations of high str.s, This Ulntultlv.ly 
correct since the error measure is related to the calculation of the 
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•trains and hanca proportional to the stresses. Finally, Figure 5 
shows the adapted aesh and the analysis results on that sash. In 
general the error Indicator has equilibrated to reasonable values and 
the answer Is considered to be uniformly reliable In the domain o£ 
Interest for the number of elements used. 

Conclusions 

Ve have presented In some detail the theory and development 
behind a three-dimensional hexshedral finite element mesh generator 
vorlclng directly from a solids model database. In conjunction with 
the mesh generator are co-developments In automatic error analysis an 
adaptive grading to produce uniformly reliable analyses. 

Research and development pertaining to the overall system Is 
continuing. Development of the mesh improvement schemes and inclusion 
of more topologically complex finite element primitives is proceeding. 
The mesh generation phase is only automatic in each finite element 
primitive. Full automation Is Impossible without automating the 
finite element primitive decomposition. Research is underway in the 
general area of feature recognition as applied to the process of 
primitive decomposition. In addition, work is continuing in the areas 
of automatic error analysis and adaptive grading. This work is 
primarily seen to remain in the area of adaptive grading because of 
the predominant nature of the commercially available finite element 
analysis codes. 




Figure 2. The Solid Model of 
the Primitive 
Decomposition 


Figure 3. The Minimum Principal 
Stress Contours on the 
Hexahedral Mesh 
Generated. 
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Figure 4. Error Contours from the 
Automatic Error Analysis 



Figure 5. The Minimum Principal 
Stress Contours After 
Adaptive Gridding 
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